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1.0 


INTRODUCTION 


Over the past several years, much effort has been expended in an 
attempt to size the power requirements for an aeronautical satellite system 
wnich would operate in the VHF oand. One reason for this effort is the 
uncertainties surrounding the variability of the signal levels due to the 
propagation path. In particular the allowances that should be made for 
absorption, scintillation and multipath fading have not been ascertained 
wi*h a high degree of confidence or agreement. Additionally, the degree to 
which these same propagation phenomena might affect ranging and position 
fixing accuracy has not been established. 

This report has been prepared to summarize the current state of 
knowledge regarding tropospheric and ionospheric scintillation and multipath 
at the VHF frequencies of interest for an aeronautical system serving the 
north Atlantic oceanic region. Only the effect of these factors on signal 
level fluctuations will be considered, and this will be done without involvement 
in aeronautical operational procedures or specific aeronautical telecommunica- 
tions systems designs. 
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2.0 


SCIOTIUATION 


2.1 


TROP0E i I • :.R IC JC LNT 1LLATI CM 


"hen satellites are viewed at very low elevation angles, the ray 
path traverses the stratifications of the earth's tropospheric refractive 
in lex profiles suer, that fading results due to the synoptic meteorological 
patterns defocusing the lobe pattern of the antenna and/or multipath prop- 
agation in the troposphere. Such fading, having fade depths of 5 to 30 db 

or i.*ore , has been experienced on :;iany line-of-sight microwave paths and is 

) 

kno\/n to be essentially frequency independent. 

VKF tests employing the »\TJ-1 satellite at an elevation angle of 

5 ' experienced higher and more frequent fade depths than might have been 

expected from purely ionospheric irregularities at a mid-latitude station 

( 7 2 ) 

m Maryland. “ There appeared to be minimal correlation of scintillation 
with s oread -F and little or no correlation with the local ..-agnetic activity 
Index. It is suggested that this data might have contained a combination of 
tropospheric and ionospheric scintillation. 

( 10 ) 

.ropospheric scintillation was observed by Tadiiio using the 
run as a signal source at four microwave frequencies: 3800, 4 995, 2595 anu 
1415 MKz. These observations are somewhat indeterrinate in evaluating the 
magnitude and time distribution of terrestrial tropospheric scintillation 
since various point sources on the solar disk vary in intensity, size and 
number producing variations in solar emission characteristics. Further, 
although 3adillo reports that frequency dependence varies in an undetermined 
manner, there is no evidence in these tests to refute the observations on 
terrestrial microwave systems that tropospheric scintillation is essentially 
Independent of frequency. The multiple source characteristics of co-ar 


I 


/ 


emission could explain the variations in frequency dependence observed. 

One significant conclusion drawn by Dadillc is that "ocintiilation decreases 
with elevation angle, becoming negligible at and greater." 

Since it is generally agreed that an aeronautical satellite 
system would not attempt to use angles below 10° elevation, and that the 
aircraft will generally be above the troposphere, it is not likely that 
this sort of signal fluctuation will be a problem. 
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2.2 IONOSPHERIC SCINTILLATION 

2.2.1 General 

The term scintillation has been applied to the amplitude and phase 

variations resulting from a redistribution of the signal amplitude caused by 

scattering of a radio wave passing through irregularities in the electron 

density distribution in the ionosphere. Absorption plays no part in 

scintillation. Amplitude variations result from the recombination of waves 

( 2 ) 

that emerge from the ionosphere with changed phase. 7 The ionospheric 
irregularities responsible for scintillation are generally much larger than 
the wavelength (especially at VHF); the usual range reported for either 

( 3 ) 

correlation or peak-to-peak shadow being from a few tenths to several kilometers. 

Ihe height of the irregularities are generally agreed to be in the 

F-region^ of the ionosphere with perhaps an additional . -region 

involvement at middle and high latitudes. 

The phenomenon of scintillation has been analyzed by use of models 

based on phase changing screen diffraction theory. while this approach 

is useful in attempting to predict or confirm some of the characteristics and 
dependencies, it has not provided a basis for practical communications system 
engineering. What is needed is an empirical approach which statistically 
describes the scintillation in terms of fade depth, duration, and frequency of 
occurrence with geographic location and other factors as parameters. Such an 
approach has only partially been started. There is considerable data on both 
radio star and satellite scintillation, but nearly all of it has been reduced 
in terms of "scintillatior index" which has a number of 

different mathematical meanings depending on the person using it. None of the 
definitions however, can be directly related to cumulative distributions desired 
by the coonnuni cat ions engineer. The information contained in the following 
paragraphs is based primarily on empirical scintillation indicies. 
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2.2.2 latitudinal and Longitudinal Dependence 

One of the dominating factors in scintillation is geomagnetic 
latitude. Aarons et al depict the ionospheric irregularity at night, 

which roughly corresponds lo the maximum scintillation regions of the world, 

( 6 ) 

as shown in Figure 1. The density of the hatched areas roughly represents 

fading depth. The high latitude region is better illustrated in Figure 2 

which is an update (to be published) of the high latitude irregularity 

structure during quiet magnetic conditions presented in reference^ K 

In order to obtain a more specific picture of the dependence of 

scintillation on geographical regions, it '/ill be necessary to resort to 

(3) 

observations made at a number of widely separated locations. Allen and 
Aarons et al^have provided a start on this by analyzing scintillation index 
at stations in the high latitude, mid-latitude end equational regions of the 
world TMs work has indicated the following: 

a) High Latitudes 

Sci it illation is present over the entire polar cap with some 

evidence of lower ratee and amplitudes over the northern polar 

( 6 ) 

cap than over the auroral oval. 

b) Mid- Latitudes 

Significantly less scintillation is present in the mid-latitudes 
between the lower ooundary of the sub-auro.-al region (about 55 
tc b0° GML) and the upper boundary of the equatorial zone 
(about 10°-15° CML^^ Near these boundaries, however, large 
fluctuations in index are possible and the ooundary edge is 
jagged, especially during magnetic disturbances. 

c ) Equatorial Latitudes 

Tne longitudinal extent of the irregularity region considered 
responsible for scintillation is a function of time of day in 
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boundary edge jacked 
during both ouiet 
period? and magnetic 
disturbances. 


'GEOMAGNETIC 

TIME 


O DYSON (1969) UPPER BOUNDARY (MEDIAN 
VALUES) Kp < 3. (alouette 2, topsice probe) 

X DYSON LOWER BOUNDARY. 


— SCINTILLATION BOUNDARY Kp r = 0,1 . 

kARCH. ,bT AL. (y UIET MAGNETO CdiDITIONS) 


Figure 2. High 'atitude Irregularity Structure 
(Courtesy J. Aarons, Private Communication, to be published 
in Journal of Cecphysical Research) 
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the equatorial regions with the uaxisdim concentrated around 
local midnight along the geomagnetic equator (See Fig. l). 

2.2.3 Diurnal Dependence 

At high latitudes, high scintillation Indices are prosent over many 
hours of the day^ with a definite diurnal period. Minimum scintillation 
occurs during local daytime with a pre-midnight maximum as shown in Figure 3. 

It may also be noted that for this data reported from Thule, an anti-correlation 
exists between fading period and fading amplitude. 

At mid-latitudes the diurnal variation favors a post-midnight 

M 

maximum but there is a lower probability of occurrence than the other regions. 

In the equatorial regions, the diurnal pattern is similar to those 

at high latitudes with a pre-midnight maximum; however, equatorial scintillation 

( 6 ) 

tends to start more abruptly reaching maximum fading in a few minutes. 

2.2.4 Seasonal Dependence 

Information on seasonal dependence is scarce. There seems to be 
general agreement that a seasonal peak occurs at the equinoxes in the equatorial 
zone ^Seasonal dependence at other latitudes is less clear. 

2.2.5 Solar Dependence 

The effect of increased solar flux on the irregularity picture is 
to increase the scintillation index on the geomagnetic equator and to decrease 
the latitudinal extent of the scintillation region A comprehensive study 
of the variation of scintillation on 38 MHz observations of Cassiopeia A at 
Cambridge, England indicated a striking correspondence between mean annual 
index of scintillation and twelve month running mean sunspot numberf^^ 
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Magnetic \ctivitv Dependence 

One of the most pronounced effects on scintillation is 

associated with disturbances of the earth's magnetic field. At high 

and mid- latitudes , there is a strong positive dependence on K indices. 

However, at the equator, there is conflicting information. Studies 

for the African continent have shown a negative correlation with higher 

(o) 

K indices in I960 while data from Peru in 1967-1969 does not corrob- 
orate this,^^ 


2.2.7 Frequency Dependence 

Frequency dependence of scintillation appears to be a function 

of several parameters. Analytically, using the phase screen diffraction 

model, the parameters include the thickness of the screen and the viewing 

distance from the screen. For a thin screen and a viewing distance much 

2 

less than the Rayleigh distance (Z^ = L where Z^ is the Rayleigh 
distance, L is the representative size of the ionospheric irregularity, 
and A is the operating wavelength; scintillation is nearly inversely related 

( 3 ) 

to the square of the frequency. '.Then the viewing distance j.s greater than 
the Rayleigh distance, the depth of scintillation becomes inversely related 

(3) 

to the first power of the frequency. Empirically^ Allen reports that for 
low depths of scintillation, the dependence was the inverse square law; but 
as the depth of scintillation increased, the dependence flattened. 
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2.2.8 llevAtion \nrle Dependence 

The variation of amplitude scintillation with elevation angle is 

( 2 ) 

not well known. Lawrence, et al predicted a dependence of the cube of 
the secant of the zenith angle employing plane-ionosphere geometry and 
assuming the irregularities to be spherical in shape. They report the 
results of several workers which indicated that the increase with zenith 
angle is appreciably less. 

( 13 ) 

Hriggs and Parkin theoretically treat the problem for both 
radio star and satellite cases. Their equations predict (for the radio 
star case) a square-roct-of-the-secant of the zenith-angle dependence when 
the observer is in the near zone, i.e. when: 



r 

o 


wavelength 

distance from the observer to the irregularity 
the distance where the correlation falls to l/e in the 
direction of the transverse field (take as 1 km for 
irregularities in the F-region) 

However, for the synchronous satellite case, the distance to the source (2^) 
must be considered. For such a case, employing a wavelength of 2 meters 
(150 MHz) the observer is not in the near zone such that the dependence 
becomes a more complex function of the zenith angle and the distances 
between the observer and the irregularity (Z^) and between the irregularity 
and the satellite t^). 


where 

a 

z. 
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To date empirical data is both sparse and divergent and one 
estimate is that the indices increase "more or less" as the secant of 

(4) 

the zenith ancle. Figure 4 shows a comparison anonf the various 

correction factors that have been considered. Currently there is insufficient 
evidence that Any of these may be used with a high degree of confidence. 


2.2.9 Magnitude and Statistics 

2. 2. 9.1 High Latitude 

Observations of ATS-3 at 137 mHz at Thule, Greenland (76.6*N GHL) 

showed fades of 3 to 15 db over one day during a magnetically quiet period 

with fade periods from 2 to 60 seconds (see Figure 3). During an intense 

magnetic storm, higher amplitude scintillations were observed. Fades of 

(14) 

10 db or more are normal daily occurrences. This data has not been 

processed in terms useful in establishing requirements as a function of 
percent time availability. Additionally, the elevation angle to the satel- 
lite was 3.7° which may have permitted tropospheric scintillation to con- 
taminate the results. 


Observations at 136 KHz on non-synchronous satellites at Fairbanks, 
Alaska showed somewhat less but still significant scintillation.^^ Figure 
shows that for three years prior to 1968, an average of 25 percent of the 
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a Elevation Angle (Degrees) 

Figure 4. Elevation Angle Correction Factors for Ionospheric Scintillation 
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Figure 5 . Amplitude Distribution of Peak-to-Peak Scintillation of 
136-MHg Satellite Signals at Fairbanks, Alaska 
(Reproduced from Coates & Golden, Ref. 15 ) 
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data acquisition passes had peak-to-peak scintillations of 12 db or greater. 
During the most disturbed month (September, 1967), 65 percent of the passes 
had scintillations greater than 12 db and 35 percent had scintillations greater 
than 30 db. 

AT5-3 VHF signal strength recordings at Narssarssuaq, Greenland 

consisting of 3581 one hour samples covering September 1968-Augu*t 1969 show 

the higher percentage of occurrence at the higher geomagnetic Latitudes 
It can be seen from Figure 6 that scintillation indices of 90% or more 

(equivalent to P^- of 12.5 db) occurs about 20% of the time during 

the worst hours of the day at this location. This should be compared to 

equivalent data taken at Sagamore Hill, Mass, which may be considered a 

mid-latitude station (see Figure 7) where the severity seems to be only 

one fourth as great. 

2. 2.9.2 Mid-Latitudes 

Many satellite and radio star scintillation observations have been 
reported at mid-latitude stations but unfortunately most of the data has been 
reduced in terms of scintillation index. For this reason, the results will not 
be discussed in detail. Information at several widely separated locations may 
be found in the following references: 


Randle Cliff, Md. 

Ref. 

(12) 

Sagamore Hill, Mass. 

Ref. 

(17) 

Greenbank, W. Va. 

Ref. 

(17) 

Boulder, Colo. 

Ref. 

(18) 

Urbana, Illinois 

Ref. 

(19) 

Baker Lake, Calif. 

Ref. 

(19) 

Stanford, Calif. 

Ref. 

(19) 
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FREQ. 137 MHz 
SEPT. 1968- AUG. 1969 
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18 24U.T. 24 03 06 C 

LOCAL MEAN TIME 

Percent Occurrence of Scintillation et *Urss«xs8u*q on AT5-3 Slgnale 
(Reproduced from Whitney, Ref. 16) 
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2.2. 9.3 Equatorial Latitudes 

The most extensive statistical scintillation data available in the 
equatorial region was taken o.» ATS-3 at 137 raH* and is .'nown in Figure 8. 

The more rapid onset of severe equatorial scintillation can readily be seen by 
comparing these curves with Figure 6 for the high latitude station at 
Narssarssuaq. Note that 90£ scintillation index occurs about 30£ of the worst 
hours of the day as compared to 20£ at Narssarssuaq. However t the gross scin- 
tillation throughout the day is actually less at Huancayo due to the more rapid 
build-up and decay. 

Further evidence of the severity of equatorial scintillation is 
shown in Figure 9 where the percent of scheduled minitrack passes missed 
during a 12-month period at three NASA stations is shownf^^ The highest 
activity is at Lima which is on the geomagnetic equator. During March 1967 
at this station 68 percent of the scheduled passes were missed during the worst 
hour between midnight and 0100 local timef^^ 

One of the significant reports of equatorial scintillation resulted 

from actual ATS-1 satellite communications support used for backup on the 

Apollo 11 recovery ship U.S.S. Hornet and the tracking ships Huntsville, 

( 21 ) 

Mercury, and Redstone. 

During the 12 day period reported, these ships were within about 

+20 degrees of the geographic equator and with the exception of the Redstone, 

all were within 30° longitude of the ATS-1 subsatellite meridian. Thus, 

relatively high elevation angles prevailed during the tests. Figure 10 shows 

a sample of the time distribution of scintillation when the Hornet was at 
o 

about 11 (about 15° invariant latitude). It depicts evenly weighted 

data over a 24 hour period and cannot be construed as the average scintillation 
experienced at that location. However, it does confirm that high scintillation 
indices may be expected in the vicinity of the geomagnetic equator a significant 
percent of time. 

18 
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Figure 9. Scheduled Minitrack Passes Missed by Hour of Day 
During 12-Month Period Due to Propagation Distortion 
(Reproduced from Coates & Golden, Raf. 15) 
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2. 2. 9.4 Cumulative Distribution of Scintillation 

The magnitude and statistics of scintillation index (SI) presented 
in the preceding paragraphs are not particularly useful to the systems design 
engineer since they do not describe the percent of time that a system may be 
degraded. 

Tables 1 thru 4 have been prepared in order to obtain a ''birds-eye" 
view of the useful data on scintillation that have been reduced to date (i.e., 
those which have been processed in terns of cumulative distribution). It can 
be seen that only two locations (Hamilton, Mass.* and Churchill, Canada) repre- 
sent data directly reduced to cumulative distributions. Table 4 shows data 

converted to cumulative distributions from tabulations of 15 minute scintil- 

(22) 

lation indices . The technique involves the construction of six model 
distributions representing six quantized levels of scintillation index. Then, 
the percent occurrence of each of these quantized levels is determined for the 
various stations, satellites and time periods for which scintillation index 
data is available. The cumulative distributions are then computed from these 
models and percent occurrences. This conversion technique has the potential 
of making available much additional data in terms of cumulative distribution 
v/ithout the need to reprocess the original raw data; however, further effort 
is required to verify some of the assumptions both explicit and implicit in 
the technique and to confirm the ability to extend the technique to stations 
other than those for which the models were constructed. 

The largest group of statistical data shown in the tables were ex- 
tracted from Intelsat I and Intelsat II. The Air Force Cambridge Research 
Laboratory (AFCTL) monitored the 136 MHz telemetry beacons of Intelsat 1 for 2280 
hours of useful data between 5 May and 9 August 1965 and of Intelsat II (F-3) 


-^-Otherwise known as Jagamore Hill, Mass. 


or 


• O • • X* 


Cambridge Research Lab. 
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for 2*050 hours of useful data between 11 May 1967 and 17 February 1968. 
These data, which were reduced by COMSAT* and published by the CCIR 
(reference 23), are shown reproduced in Figure 11. The data were taken 
at Sagamore Hill, Mass, (about 53° invariant latitude) with Intelsat I 
at an eleyation angle of 25° and Intelsat II at an elevation angle of 
13 

^COMSAT has been reverifying their data processing techniques. There is 
indication that the curves may be changed considerably. 
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TABLE 1 - INTELSAT DATA (CCIR REPORT IV/IO 67 ) 


Station 

USAF Cambridge Research Lab, 

, Hamilton, Moss. 

Geog. Lat. 

42*47* N 

Geog. Long. 

70*51' W 

Geom. Lat. 

53°N 

Satellite 

INTELSAT -I 
EARLY BIRD 

INTELSAT -I I 
CANARY BIRD 

Sat. Long. 



Freq. 

136 MHz 

136 MHz 

Elevation 

e 

ITV 

CVI 

13* 

Date of Data 

5 - 5-65 to 8 - 9-65 

5-67 to 2-63 

Processed Data 
Hours 

2280 

4050 

Only SCIN 
Data Reduced 



All Data Reduced 

X 

X 

Sunspot # 

15 (See Note *A) 

33 (See Note *A) 

Magnetic Act. 

6 (See Note *B) 


95* 

-0.2 dB 

-1.2 dB 

9 9* 

-0.8 dB 

-4 .8 dB 

99.9* 

-3.0 dB 

*c 

99.99* 

*C 

*c 

NOTES 

1. Azimuth 126.7° 

2. Sub -ionospheric coordi- 
nates 38°N lat., 64°W long. 

3. Column 4 composite of columns 
1 , 2 & 3 plus non scintil- 
lation data 

4. Scintillation data less than 
2* of all data recorded 

5. Data reduced by COMSAT 

1. Scintillation data, 17* 
of total recorded data 

2. Data reduced by CCMSA? 

References 

24 

34 


- Sunspot number Is average number of sunspots over period of data collection 
(yearly running sunspot number per CCIR Report 246-1) 

- A p number mean average value of monthly readings over period of data collection 

- Inconclusive data 
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* IN ... .-u> DATA (l .... 


3'. ati n 

’3*‘. :ambrl lge . ( Search lu»b., ilamiit. n, Mass. 

DEI 

42 +7* :j 

Geog. * ng 

70"5i'v 

IS 

53 e N 

3a belli :e 

▼ .n CAI' t 

- *» 1 .4 .Oil 4 "* 

iA . ' 

-11 

A.i/i • i D 

3a t . L ng . 



a- 

(U 

l3^> MHz 

136 a » 

Levati >n 

25 

13 ' 

Date of 

Data 

5 - 5-05 i 
to 

6 - 10-65 

to 

7-17-65 

1 

7-19-65 

to 

■ - } -65 

5 - 11-07 

to 

6 - 10-6 7 

o-lo -07 
7-10-67 

7-10-67 

to 

'- 10-67 

1-10-67 

9-L0-67 

9 - 10-07 10.10-67 
to to 

10 -10-6/ i.1 -10-oY 

i ned 

Data :! irs 

21.2 

.7 


156 

235 

205 

. . 

• 4 . 0 


only 3CIj.' 
Data .\educt d 

A 

X 

X 

X 

X 

X 

X 



All Data 
Seduced 








X 


S IQfl ' 

See Note */ 

15 

15 

16 

n 

67 

92 

107 

77 


Magnetic 

Note *1 


10 

8 

25 

12 

CJ 

O 

) 

16 

10 

>5 1 

1 


-1.5 U 

-3 • • 

. 

. • : 

-l.i cLl 


- . • . . 

99 % 

- . 

BSE 

-3.6 cLB 

. 

-3.5 d.: 

-4 , 4 d»S 

. 

H 

. : 


- .. 

-7.3 cd 

-7.5 dB 



- . 

- . 



>9.99 

- .< : 

* ' 

*■(' 

*c 


1 

. 

N * 

-10 db 

-0 . 1 LI- 

IJOT3S 

1. Azimuth 126.7 

2. 3 .b -ionospheric 

coordinates 3 ’1 

at. 64 W long. 

3. Column 4 com posit 
of columns 1, 2 & 

3 plus non -scintil- 
lation data 

4. Scintillation data 
less than 2}> of all 
data recorded 

5 . Data reduced by 
CLMaAl 

ditl| 

2. Data Reduced by COMSAT 

References 

24 

34 


- Sunspot number used from month with most nunber of data collection lays (yearly 
running sunspot number per CCIR Report 246-1) 

5 - A 3 number, mean value of daily reading for month with most :. miber of oata collect! n 


•r - 

- Inconclusive data 
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TABLE 3 - LES 5 & 6 DATA 
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Station 

Churchill 

Geog. Lat. 

58.8 *N 

Geog. Long. 

94.2°w 

Geom. Lot. 

68.4* N 

Satellite 

LES 6 

LES 5 

Sat . Long 


Drifter 

F req . 

254 MBs 

228 MHz 

Elevation 

15 ° to 22 ° 

0 ° to 26 ° 
17° Mean 

Date of Data 

8 - 28-69 to 
10-4-69 

8 - 28-69 to 

10-4-69 

Processed Data Hours 

741 

115.2 

Only SCIH Data Reduced 



All Data Reduced 

X 

X 

Sunspot # (See note *A) 

106 

106 

Magnetic Act. 



) p i (See note 1) 

- 3.6 dB 

- 6.8 dB 

(See note 1) 

-9.5 dB 

-12 dB 

99.90 (See note l) 

-20.4 dB 

-22 dB 

99.99/j (See note 1) 

-24 dB 

-27 dB 

NOTES 

1. Data scaled to 136 MHz by l/f 

2. 98 1> of data recorded showed SCIN 
activity 

References 
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-* 

A - Sunspot number used from month with most number of data collection 
^ys (yearly running sunspot number per CCIR Report 246-1) 
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!'AHLE 4 - INTELSAT fit ATS DATA (REDUCED BY AFCRL) 


Station 

CAKBRIXE 
RESEARCH LAB. 

KARSSAItSSUAW 

GREENLAND 

HUANCAYO 

PERU 

Geog. Lat. 

42 47’;: 

60 °N 

11 i 

Oeog. Lat. 

70*51*W 

45 W 

rrv 

Georn. Lat. 

53 • 

7i.5°u 

0° 

Satellite 

INTELSAT IT 

ATS -3 

mM/M 

ATS -I 

Sat. Long 


45°W 

45°w 

150 U W 

Freq. 

136 rdHz 

137MHz, 

137MHz 

136 Mz 

elevation 

iy 

• 

20 " 


Date of Data 

5 - 11-67 to 
2 - 17-66 

7 • 

9 - 29-69 

9-17-63 to 
8-14-69 

7-14-6? to 
2-2 -69 

Processed Data 

Hours 

4212 

12330 

3314 

6513 

>nly 3c in Data 

Reduced 





All Data Reduced 

X 

X 

X 

X 

s inspot f (s ee 

Note # 

A) 3 

•107 

lOo 

104 

Magnetic Act. 





95 % (See Note *B) 

1 -0.6 dB 

* -0.4 dB 

-2.2 dB 

-1.5 dB 

n *» 

22k 

ti 

-1.6 dB 

-1.4 dB 

'•5.3 dB 

-5.3 dB 

tt w 

99.9-i 

• 

-4.6 dB 

-5.0 dB 

-12 dB 

-12 dB 

, it fi 

>9.99* 

it 

-10.0 dB 

-10.0 dB 

-K- - 

* • 

notes 


1 . l'o848 
15 -minute 
SCIN values 

1. 49421 
15 -minute 
SCI II 
values 

1 . Dirdng above 
time ATS -3 
m^vea from 

94 W to 45 °W 

1. 2oo?4 
15 -min lie 
SCIN 
va lues 



2. Data re- 
duced by 
AFCRL 

2. Data re- 
duced by 
AFCRL 

2. 14255 

15 -minute 

SCIN 

values 

2. Data re- 
duced by 
AFCRL 





3. Data re- 
duced by 
AFCRL 
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* 

A - Sunspot Number is average number of sunspots over period of data collection 
(yearly running sunspot nimber per CCIR Report 246-1) 

*B - Converted from statistics on occurrence of scintillation index (conversion 
technique not yet verified) 

* t 

C - Inconclusive data 

































































The 1965 data showed that 90% of the fades exceeding 3 db lasted 

less than 25 seconds while the estimated median duration of 6 db fades was 

( 2 L) 

about 5 seconds. In the 1967 data, 5*# of the 3 db fades and 63% of 

the 6 db fades lasted for less than 10 seconds. 

A significant difference can be seen in Figure 11 between the 1965 
and 1967 data at the ’’tails" of the distribution curve. The reason for this 
has not been confirmed but one can speculate on some possibilities. First, 
assuming the elevation angle dependence to vary as the secant of the zenith 
angle, the difference between the two curves might be expected to be slightly 
more than two to one (in db). The observed difference, even at 99% t is far 
greater than this. Second, the sunspot number was in the order of 15 for the 
period the data was collected in 1965 and in the order of 83 for the period the 
data was collected in 1967-1968. According to Aarons, ^ the logarithm of the 
depth of scintillation depends nearly linearly on the yearly average of the 
sunspot number. Thus, it is more reasonable to account for the differences in 
this manner. However, the divergence of the curves implies that even both these 
dependencies cannot fully explain the difference. Further, it has been noted 
that on May 25 and 26, 1967, a severe magnetic storm occurred. When the 
distribution curves for individual months was investigated, it was found that 
more than half the deep fading responsible for the rapid fall off of the 1967 
data occurred in the month of May 1967. This may have also contributed to the 
differences noted above.* 

For convenience in analysis, the data in Figure 11 have been tabulated 
in terras of four time periods 95, 99, 99.9 and 99.99% of the time in the tables. 

In addition, the monthly scintillated data that made up these curves have been 
tabulated. Although some trends with regard to geometric dependence, solar 
and magnetic activity may be observed, there is no clearly deterministic relation- 
ship evident. The Intelsat data taken at Hamilton indicate that 
♦Data processing methods may also have contributed to the differences. 
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about 5 db fading for 1# of the time may be experienced at 13° elevation 
angle in the mid-latitudes about 53° GML during nigh sunspot periods and 
about one db may be experienced for less than 1? of the time at 25° ele- 
vation angle during low sunspot periods. 

The Churchill data taken or. LES-5^*^ and LI'S-6 ^ shown in Table 3 was 
scaled from 228 and 254 MHz, respectively to 136 MHz by the inverse frequency 
law. Thus, the data presented in this table may be somewhat suspect when com- 
pared to data taken directly on 136 NHz. However, it does agree with the 61 
observations that greater fading occurs at the higher geomagnetic latitudes. 

Considering that Tables 1 thru 4 represent all the data known to the 
author to date that ha* been reduced to cumulative distributions, it is evident 
that we are severely lacking in ionospheric scintillation data in a form use- 
ful for determining systems design parameters, especially in the high geo- 
magnetic and equatorial latitude regions. 

2.2.10 Geo grap h! c Ti . stri buti on of Scintillatio n 

A more or less qualitative picture of the geographic distribution of 

scir illation fading was summarized earlier in this paper. At present this 

seems to be the most significant geographic distribution available. 

Many concepts have been considered for extrapolating the little 

(3) 

known data to other geographic locations; however, Allen'^ 'cautions against 

attempting to use empirical relationships such as geometric dependence, in regions 

(13) 

other than those from which they are derived. More recently, Aarons has 
suggested that one might consider scintillation to be relatively constant along 
lines of invariant latitude in the high (and perhaps) mid- latitudes. However, 
there does not appear to be a sufficient amount of supporting data to apply 
this theory except on & very "gross" basis. 
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Since a truly representative geor-aphical picture of scintillation 

fading is not available, a general idea ol its impact on a VHF aeronautical 

system may be obtained from the experience of actual airlines tests during a 

major magnetic storm. During vbe period of October 24-31, 1968, airline 

tests were conducted employing ground stations at Annapolis, Md., Chicago, 

111., and Wichita, Kansas in conmunication with three flights between the 

(26) 

U.S. and Europe through the ATS-3 satellite. 

During this period, which encompassed a major solar flare and associ- 
ated geomagnetic activity (Max. Kpp * 8), 24 db peak-to-peak scintillation was 
experienced at the Annapolis and Chicago ground stations. However, the Wichita 
station displayed none of the variations reported at the other stations. 

The aircraft scintillation data was not analyzed statistically since 
it would have reflected many other propagation factors, however, it was reported 

that the link was available 8251 of the time and during this time usable voice 

( 26 ) 

conmunication was achieved 59# of the time. The highest latitude reached 

(on one flight) was about 62° invariant latitude and most of the flight courses 
were above 55° invariant latitude. Thus, the aircraft was experiencing a mix 
of mid-latitude and "sub-auroral" scintillation. These aircraft tests sh<w that 
some consideration must be given to scintillation effects in a practical operating 
case. 


(27 ) 

An empirical model has been proposed by Fremouw and Bates v ' consisting 
of three additive terms corresponding to equatorial, mid- latitude, and high 
latitude irregularities. Such a model could provide the necessary information 
for the systems designer provided a sufficient amount of empirical data cjuld 
ha made available to establish the coefficients within the equations and the 
results presented in terms of fade statistics useful to the engineer. 
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2.2.11 Discussion 
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Although there is a large amount of data available on scintillation, 
meaningful statistics in a form useful to the system engineer is available 
from only one mid-latitude location (Sagamore Hill, Mass.). Some limited 
VHF data in useful form is slowly becoming available from other locations, 
especially Churchill, Canada. Because of quantity limitations and the 
critical sensitivity of this type fading to the parameters of the ionospheric 
irregularities, extension of existing data to other condiLion3 (such as dif- 
ferent elevation angles, other geographic locations, phases of the sunspot 
cycle, higher frequencies, etc.) is questionable. On a gross basis, the best 
estimates at this time would suggest the following scintillation index depen- 
dencies: 

1/2 

Elevation angle - (secant of the zenith angle) 

Frequency - inverse frequency 

Sunspot number - logarithmically with yearly average sunspot number 
Geographic - constant along given geomagnetic latitude 

It is quite difficult at this time to draw any firm conclusions 
regarding the statistical distribution of ionospheric scintillation or to 
extend the available data into geographic areas beyond those in which it was 
taken. Estimates of margin allowance for scintillation at this time can only be 
made from the limited data contained in Tables 1 thru 4 of paragraph 2. 2. 9.4 end 
they would primarily apply to a limited area in the mid-latitudes. Since judg- 
ments regarding aeronautical routing, operational procedures and aeronautical 
tele-comnunications systems design must be considered in interpreting this data, 
it is beyond the scope of this paper to suggest a specific system margin require- 
ment. 
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POLARIZATION AND FARADAY ROTATION 
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When a linearly polarized wave penetrates the ionosphere, its vector 

angle in space is rotated by the interaction between the wave and the geomag- 

( 2 ) 

netic field influenced by the integrated electron content along the path. 

This (Faraday) rotation has a strong diurnal characteristic as a result of the 
diurn variation in total electronic content. At VHF frequencies, the difference 
in rotation between the daytime and nighttime ionosphere may amount to several 
revolutions and, thus, must be considered in an aeronautical satellite system 
operation at these frequencies. 

As a result of the changes in the geometry of the ray path through 
the earth’s magnetic field and the variation of the field intensity at different 
latitudes, the rotation will exhibit a latitudinal dependence. Such a depen- 
dence was reported in measurements employing ATS-3 on the German research-ship 

(28) 

"Meteor" where the number of daily rotations increased as the ship sailed 
north from 14.5° to 21.5°. 

When linearly polarized antennas are employed at both the satellite 
and the aircraft, severe fading amounting to 30 db or more may result when 
the polarization of the received signal is orthogonal to the polarization of 
the receiving antenna. This problem can nearly be eliminated if circularly 
polarized antennas are employed at bo*h ends of the link; or even if only one 
antenna is circularly polarized, if one is willing to accept a fixed 3 db signal- 
loss. 

The degree to which the effects of Faraday rotation may be reduced 
depends on the achievable axial ratio (or circularity) of the polarization of 
the antenna designs. An axial ratio of 2 db in a VHF antenna is reasonably 
achievable while the VHF aircraft antenna might vary between Z and 10 db 
depending on the design and "look angle" to the spacecraft. 
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The axial ratio of the aircraft antenna is also a parameter in the 

analysis of multipath fading due to the variation in reflection coefficient 

as a function of the polarization angle of the incident reflected wave. 

(29) 

Levatich has analyzed this problem and presented curves which will be 

presented and discussed in the "Multipath" section of this paper. 

It is generally agreed that VHF polarization loss (neglecting 

multipath) would generally be between 1 and 2 db* assuming a circularly 

polarized spacecraft antenna and a reasonably well designed (elliptically 

polarized) aircraft antenna. 

4.0 EARTH REFLECTED MULTIPATH 

In satellite-to-aircraft conmuni cat ions links, amplitude fluctuation 

of the received signal due to addition of out-of -phase signal components reflected 

from the earth to the direct signal mr \j cause serious rapid frding. The magnitude 

->f the fade depth is dependent on a complex set of variables , the predominant 

ones being: 

a.) The aircraft antenna discrimination between direct and 
reflected signals. 

b) The reflection coefficient of the reflecting surface. 

(This makes sea state an indirect variable. ) 

c) The elevation angle at which the aircraft views the satellite. 

d) The ellipticity ratios of the satellite antenna and the 

aircraft antenna in the direction of the different propagation paths. 

e) The polarization of the incident wave on the reflecting 
surface. (This makes Faraday rotation by the ionc sphere an 
indirect variable. ) 

f) 't’he altitude of the aircraft. 

♦depending upon aircraft attitude 
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In addition there are several variables which affect the rate of 
fading, the more significant one being the aircraft velocity. The velocity 
of v:aves on the surface of the sea at the reflection point is small enough 
to be neglected in most cases. 

The multipath problem has been examined theoretically by several 

( 29) (30 (31) (32) (33) (29) „ , 

authors. Levatich presents curves of absolute 

depth of fading as a function of elevation angle with aircraft antenna axle' 
ritio and three values of aircraft antenna discrimination as parameters. 

Figure 12 shows these curves where it can be seen that at elevation angles 
of 10° and an aircraft antenna discrimination of 3 db, about 7 db fade depths 
may be expected, 3ince these curves assume a two-ray case and smooth sea 
conditions, they probably approximate a realistic worst case fade depth. 

xperimentally, a series of tests were conducted with the VHF trans- 
ponders on the AT3-1 and IT 1-3 satellites employing a calibrated test spectrum. 
One of the results of these tests was a three dimensional fading model with 
amplitude, time and frequency as coordinates. This model demonstrated a 
selective fading pattern within the 30 KHz test spectrum. It also showed the 
fade periods to be in the order of 5 seconds for the particular set of flight 
conditions. 


(33) 


Fro:, a practical systems engineering point of view, the probability 
distribution of fading obtained in these tests may serve as an empirical 
guide to the fade margin reouirements for a VHF system. Figure 13 shows a 
probability distribution of signal fading where the data was normalized to 
the same median. This in effect removes the longer term fading associated 
with the aircraft orientation and aircraft antenna pattern. Thus, the 
distributions may be considered to indicate primarily the statistics of 
multipath fading plus small contributions from Faraday rotation and scin- 
tillation. 
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Figure IJ. Probability Distribution of Signal Pading 
(Reproduced from Bergaann 4 Kucera, Ref. 33) 
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OTJ'ARY AMD COHCI/JS IONS 

\1 though considerable data is available on ionospheric scintil- 
lation in terms of scintillation index, reliable data taken at 13* MHz in 
terms of cumulative distributions is available from only one mid-latitude 
station. This data shovre a spread of from 1 to 5 db fading for 1? of the 
time depending on the elevation angle, sunspot number and other parameters. 
An attempt to convert data that has been reduced in ter:ns of scintillation 
index to cumulative distributions requires further verification. Further, 
the ability to extend tne limited cumulative distribution data to geo- 
graphic areas beyond those in which it was taken is highly questionable 
especially with respect to frequency scaling and elevation angle scaling. 
r t is suggested that considerable effort remains in both converting prior 
available data to cumulative distributions and acquiring new data (especiall 
in the high latitu e and equatorial regions) before the confidence in a 
system margin allowance for ionospheric scintillation can be significantly 
raised. 

Hince polarization and Faraday rotation fading is so intimately 
related to multipath fading and antenna characteristics, it seems more 
reasonable to assess this type of fading along with multipath rather than 
-ttempt to ascertain tv;o (db) additive margins. Theoretical studies have 
indicated that the absolute depth of fade may be in the order of 7 db at 
an elevation angle of 10 0 assuming a spacecraft antenna ellipticity ratio 
of 2 db and an aircraft discriminator of 3 db and axial ratio of 10 db. 
\ctual measurements on the '75 program employing a linear satellite antenna 
and aircraft antennas representative of early designs indicated that about 
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4 db or more fading was experienced l£ of the time between 30° 
and 50° elevation angles and about 7 db or more was experienced l£ 
of the time between 0° and 30° elevation angles. The data was taken 
employing ^arly aircraft antenna designs and a linearly polarized 
spacecraft antenna, therefore any conclusions drawn from them are 
probably pessimistic when considered in light of what might be achieved 
on a future satellite program. 

It should be pointed out that in establishing systems margins 
based on the propagation statistics summarized in the report, the fade 
depths should not be added db-wise but rather the statistics of each 
type of fading shouxd be combined on a time frame basis of an hour, a 
flight duration, a day, a month, a year or a sunspot cycle, whichever 
fits the operational needs of the user. 
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